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Problem  Statement 


The  basic  objective  of  this  research  project  was  to  examine  dynamic  fracture 
criteria  to  determine  if  the  material  constants  of  the  failure  process  can  be 
obtained  by  using  a  relatively  simple  experimental  configuration;  with  the  re¬ 
quirement  that  the  dynamic  failure  process  can  be  predicted  with  sufficient 
generality,  to  allow  its  application  to  impact  conditions  with  various  geo¬ 
metries  and  stress  fields. 


Study  Sequence 

To  achieve  the  above  objectives  the  stu<fy  was  conducted  in  essentially 
five  phases.  Each  phase  was  undertaken  to  examine  a  specific  aspect  of  the 
failure  process,  and  each  one  was  a  logical  extension  of  the  one  before.  This 
study  sequence  can  be  sumiarlzed  as  follows: 

(1)  Existing  failure  criteria  were  examined  to  determine  their  ability 
to  predict  failure  under  different  impact  conditions  and  geometries. 

(2)  A  combined  theoretical  and  experimental  study  was  undertaken  to 
determine  a  failure  model  appropriate  to  the  prediction  of  static 
and  dynamic  failure  in  thin  rings  and  cylinders. 

(3)  The  model  developed  in  (2)  was  applied  to  the  prediction  of  failure 
in  thick  plates,  as  represented  by  the  uniaxial  strain  plate  impact 
experiments. 

(4)  An  improved  failure  model  was  introduced  which  incorporated  contin¬ 
uous  void  nucleation  into  the  model  developed  in  (2). 

(5)  Improved  experimental  techniques  were  developed  for  thin  ring  failure 
analysis,  to  allow  the  failure  properties  of  high  strength  steels  to 
be  determined. 
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Results 


The  main  results  obtained  during  the  course  of  the  study  can  be  summarized 
in  terms  of  each  phase,  and  these  were  as  follows: 

(1)  It  had  been  the  intention  to  use  one  of  two  existing  dynamic  failure 
criteria,  and  adapt  them  for  application  to  the  thin  ring  configur¬ 
ation.  It  was  determined  that  neither  the  accumulative  damage  model 
of  the  form: 


f  2  x 

1  (0-oo)Adt  =  K 


a  >  a 


0 


or  the  void  growth  model 


exp 


an0 


could  be  used  in  both  the  plate  impact  and  ring  configurations, 
mainly  because  the  material  constants  (oQ,  X  etc.)  appropriate  to 
one  configuration  did  not  apply  to  the  other. 

Further,  it  was  found  that  the  quasi -static  void  growth  model  of 
Rice  and  Tracey  [2]  did  not  give  consistent  results  when  applied  to 
situations  where  high  strain-rate  conditions  prevailed.  However, 
as  this  model  was  not  develuped  for  dynamic  problems  these  results 
were  not  too  surprising. 
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(2)  Recent  developments  In  the  use  of  long  wave-length  plastic  stability 
criteria,  as  outlined  by  Hutchinson  and  Neale  [3],  together  with  the 
use  of  the  Gurson  yield  function  [4],  suggested  an  entirely  different 
approach  to  the  problem  of  dynamic  ductile  failure.  This  yield 
function  has  the  form 


+  2  f  •  cosh(^)  -  f2  -  1  =  o 
m 


(a) 


where  f  is  the  void  volume  fraction,  and  Y„  is  the  flow  stress  of  the 

m 

matrix  material.  The  compressibility  of  the  plastic  material  is 
introduced  through  the  relation 

f  -  0  -  f)  (e,p  +  e/  +  e3P)  (b) 

The  theory  appropriate  to  the  study  of  the  static  and  dynamic  failure 
of  thin  rings  was  developed,  together  with  a  companion  set  of 
experiments.  This  approach  gave  very  encouraging  results,  as  shown 
in  Fig.  1.  The  material  property  associated  with  failure,  the  critical 
void  volume  fraction,  was  obtained.  It  was  also  noted  that  the 
dynamic  results  were  not  particularly  sensitive  to  the  strain-rate 
sensitivity  of  the  material.  Details  of  this  work  have  been  reported 
in  [5],  [6],  and  [7]. 
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(3)  The  above  theory  when  applied  to  a  computer  program  which  had  been 
previously  developed  to  analyze  the  spallation  process  associated 
with  plate  impact,  produced  the  results  shown  in  Fig.  2.  The 
predicted  spall  threshold  obtained,  using  the  ring-generated 
critical  void  volune  fraction,  was  not  unreasonable,  but  certainly 
higher  than  one  would  expect  of  a  model  that  did  not  include  void 
coalescence.  The  scatter  in  the  experimental  data,  due  mainly  to 
the  different  methods  used  to  define  spallation  posed  a  problem.  It 
was  felt  that  the  Sandia  data  presented  by  Butcher  and  Young  [8] 
not  only  reflected  a  consistent  definition  of  incipient  spallation, 
but  also  was  presented  with  sufficient  detail  that  it  should  be 
used  as  a  reference  for  the  theories  developed  in  this  project. 

When  the  calculated  value  of  the  peak  normal  stress  threshold  are 
compared  with  the  Sandia  values  the  correlation  is  much  better,  as 
shown  in  Fig.  3.  The  better  agreement  obtained  in  this  latter  case 
is  an  indication  of  the  Importance  of  strain-rate  sensitivity  in 
the  constitutive  model.  Another  shortcoming  of  the  model  is  that 

it  does  not  reflect  the  increased  threshold  that  occurs  as  the  pulse 
duration  is  shortened. 

(4)  The  above  model  assumes  that  all  the  void  nucleation  occurs  at  the 
onset  of  yielding.  Realistically  one  would  expect  that  nucleation 
would  occur  during  the  plastic  deformation,  and  depend  on  the  history 
of  that  deformation.  The  model  was  changed  to  Include  continuous 
void  nucleation  by  altering  equation  (b)  as  follows: 
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f  =  (1  -  f)  (e-]P  +  e2P  +  e3P  )  +  F  •  eP 


where  F  is  a  distribution  function  that  depends  on  the  history  of 
the  plastic  strain.  The  constants  appropriate  to  F  were  obtained 
from  the  history  of  the  effective  plastic  strain.  The  Sandia 
experiment  which  produced  spallation  at  the  lowest  stress  value 
was  used  in  this  case  (i.e.  experiment  4  in  Fig.  4).  The  results 
obtained  from  this  improved  model  are  shown  in  Fig.  3  and  5,  and 
the  hoped-for  change  seems  to  have  been  obtained.  To  get  some 
idea  of  the  generality  of  this  model  it  was  then  applied  to 
spallation  in  a  dynamically  loaded  thick-walled  cylinder,  and  the 
results  are  shown  in  Fig.  6.  In  this  case  the  results  were  very 
good  for  one  pulse,  but  slightly  out  in  the  case  of  a  very  short 
high  intensity  pulse.  However,  it  is  possible  that  the  loading 
pulse  was  not  as  precisely  defined  in  the  latter  case.  Although  F 
was  chosen  to  be  compatible  with  the  righ  configuration  the  ring 
programs  have  not  as  yet  been  altered  to  reflect  the  change. 

(5)  The  complexities  associated  with  determining  material  properties 
under  dynamic  loading  almost  excludes  the  possibility  of  relying  on 
one  experimental  configuration  in  this  type  of  work.  In  the  case 
of  failure  properties  the  results  of  this  study  indicate  the  ring 
configuration  can  provide  valuable  data.  The  experimental  techniques 
used  in  this  study,  as  described  in  [5],  although  ideally  suited  to 
some  materials,  may  not  provide  either  the  high  strain-rates  or  the 
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loading  levels  required  to  make  this  technique  applicable  to  high 
strength  materials.  The  relative  simplicity  of  an  exploding-wire 
system,  and  the  fact  that  the  loading  pressures  are  not  required 
allowed  the  system  to  be  adapted  to  meet  these  additional 
requirements,  and  the  configuration  shown  in  Fig.  7  was  developed. 

The  appropriate  shape  of  the  water  filled  cavity  (not  shown  in  Fig.  7) 
required  to  produce  the  appropriate  uniaxial  hoop  expansion  is  the 
key  to  this  system.  This  technique  was  developed  to  the  point  of 
showing  its  feasibility,  but  no  attempt  has  yet  been  made  to  obtain 
fracture  data. 


Computer  Programs 

The  following  computer  programs  were  used  in  the  analysis  related  to  this 
project: 


(1)  "VOIDS"  :  -  a  Runge-Kutta  solution  of  the  quasi-statically 

deformed  thin  ring;  used  to  obtain  the  critical 
void  volume  fraction. 

(2)  "DYNVOID"  :  -  similar  to  "VOIDS"  except  that  it  includes  inertia 

terms  required  in  the  analysis  of  impact  loaded 
thin  rings. 

(3)  "HOSLAP"  :  -  a  method  of  characteristic  solution  which  determines 

the  stresses,  displacements,  strain  and  void  growth 
levels  resulting  from  plate  impact. 

(4)  "CYCLOPS"  :  -  a  cylindrical  version  of  "HOSLAP"  which  determines 

the  response  of  a  thick-walled  cylinder  to  an 
internal  pressure  pulse. 
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VOID  VOLUME  FRACTION,  f  =  f/f^  VOID  VOLUME  FRACTION, 


o 


UNIFORM  STRAIN,  £* 

o 


Fig.  1  Critical  Void  Volume  Criterion  Applied  to 

Thin  Rings  of  Aluminum  and  Copper,  (n  =  0.03, 
Af  =  0.01,  Cq  =  8  x  103/sec) . 
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Normal  Tensile  Stress  at  Incipient  Spall. 
(  •  Sandia  experimental  data,  ref.  8; 

-  simple  void  growth  model; 

— -  nucleation/void  growth  model) 
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Comparison  of  Spall  Threshold  Theory  Prediction 

( - simple  void  growth; 

-  nucleation/void  growth  model) 
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